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Abstract 

The dynamical model of pion electroproduction developed in Physical Review C63, 055201 (2001) 
has been extended to investigate the weak pion production reactions. With the Conserved Vector 
Current(CVC) hypothesis, the weak vector currents are constructed from electromagnetic currents 
by isospin rotations. Guided by the effective chiral lagrangian method and using the unitary 
transformation method developed previously, the weak axial vector currents for tt production are 
constructed with no adjustable parameters. In particular, the N-A transitions at Q 2 — > are 
calculated from the constituent quark model and their Q 2 -dependence is assumed to be identical 
to that determined in the study of pion electroproduction. The main feature of our approach is 
to renormalize these bare N-A form factors with the dynamical pion cloud effects originating from 
the non-resonant n production mechanisms. The predicted cross sections of neutrino-induced pion 
production reactions, N^v^, fj,~w)N, are in good agreement with the existing data. We show that 
the renormalized(dressed) axial N-A form factor contains large dynamical pion cloud effects and 
this renormalization effects are crucial in getting agreement with the data. We conclude that the 
N-A transitions predicted by the constituent quark model are consistent with the existing neutrino 
induced pion production data in the A region, contrary to the previous observations. This is 
consistent with our previous findings in the study of pion electroproduction reactions. However, 
more extensive and precise data of neutrino induced pion production reactions are needed to further 
test our model and to pin down the Q 2 -dependence of the axial vector N-A transition form factor. 

PACS numbers: 13.60.-r, 14.20.Gk, 24.10.-i, 25.30.Pt 
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I. INTRODUCTION 



It is now well recognized that an important challenge in nuclear research is to understand 
the hadron structure within Quantum Chromodynamics(QCD). One of the important in- 
formation for pursuing this research is the electromagnetic and weak N-A transition form 
factors. Such information are important for testing current hadron models and perhaps also 
Lattice QCD calculations in the near future. In particular, we can address the question 
about whether N and A are deformed. If they are deformed, is it due to the gluon inter- 
actions between quarks or due to the pion cloud which is resulted from the spontaneously 
breaking of the chiral symmetry of QCD? 

As a step in this direction, we have developed^, in recent years a dynamical model for 
investigating the pion photoproduction and electroproduction reactions at energy near the A 
resonance. The dynamical approach, which has also been pursued by several groups El 4, H, 
EL 0, El j is different from the approaches based on either dispersion relations 11, IS , 13 1 
or K- matrix method 0, 0, 0, 0] in interpreting the data. In Refs.l and 2, we not only 
have extracted the electromagnetic N-A transition form factors from the data, but also have 
provided an interpretation of the extracted parameters in terms of hadron model calculations. 
In particular, we have shown that including the dynamical pion cloud effect in theoretical 
analyses can resolve a long-standing puzzle that the N-A Ml transition form factor G^ A (0) 
predicted by the constituent quark model is about 40 % lower than the empirical value. 
Furthermore, the dynamical pion cloud is found to play an important role in determining 
the Q 2 -dependence of G^{Q 2 ) which drops faster than the proton form factor G P {Q 2 ) as 
Q 2 increases. The predicted electric E2 and Coulomb C2 N-A transition form factors, which 
exhibit ve ry p ronounced meson cloud effects at low Q 2 , have stimulated some experimental 
initiatives Similar results were also obtained in Ref . fTol] . 

Since weak currents are closely related to electromagnetic currents within the Standard 
Model, it is straightforward to extend our dynamical approach to study weak pion production 
reactions. In this paper, we report the progress we have made in this direction. Our objective 
is to develop a framework for extracting the axial vector N-A form factors from the data of 
neutrino-induced pion production reactions. In particular, we would like to explore whether 
a dynamical approach , as developed in our study of pion electroproduction, can also resolve a 
similar problem that the axial N-A transition strength calculated ljl, 2(| from the constituent 



quark model is about 30 % lower than what was extracted [21| from the data. We will also 



make predictions for future experimental tests which could be conducted at Fermi Lab, KEK 
and Japan Hadron Facility in the near future. 

To introduce our model, it is necessary to first briefly review the dynamical approach 
developed in Ref. 0,0] (called SL model). Our starting point is an interaction Hamiltonian 
Hj = Tmb^b' which describes the absorption and emission of mesons(M) by baryons(fi). 
In the SL model, such a Hamiltonian is obtained from phenomenological Lagrangians. In 
fact, the approach is more general and this Hamiltonian can be defined in terms of a hadron 
model, as attempted, for example, in Ref. 2^], or by using the vertex functions predicted 



by Lattice QCD calculations. 

It is a non-trivial many-body problem to calculate ttN scattering and 'yN — > ttN reac- 
tion amplitudes from the interaction Hamiltonian Hj = Tmb^b 1 - To obtain a manageable 
reaction model, a unitary transformation method 0, H3| is used up to second order in Hi to 
derive an effective Hamiltonian. The essential idea of the employed unitary transformation 
method is to eliminate the unphysical vertex interactions MB — > B' with tum + fnB < f^B 1 
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from the Hamiltonian and absorb their effects into MB — > M'B' two-body interactions. For 
the pion production processes in the A region, the resulting effective Hamiltonian is defined 
in a subspace spanned by the 7rN, ^yN and A states and has the following form 

Heff = Ho + V tN + v yn + ^wN^A + r 7 Af^A, (1) 

where v n N is a non-resonant ttN potential, and v 17T describes the non-resonant 'jN <-» ttN 
transition. The A excitation is described by the vertex interactions r 7 jv<-+A for the 7iV <-» A 
transition and T^n^a for the ttN <-» A transition. The non-resonant t> 77r consists of the 
usual pseudo-vector Born terms, p and to exchanges, and the crossed A term. 

From the effective Hamiltonian Eq. it is straightforward to derive a set of coupled 
equations for ttN and 7iV reactions. The resulting pion photoproduction amplitude can be 
written as 

T^{E) - + E _ m ^_^ {E) • (2) 

The non-resonant amplitude t 77r is calculated from t> 77r by 

t J7r (E) = f 77r + t 7TN (E)G 7TN (E)v~ f7T , (3) 

where Gvtv is the tcN free propagator, and t n ^ is calculated from the non-resonant ttN 
interaction v n N. 

The dressed vertices in Eq. (J2J) are defined by 

TjN^a(E) = r 7 jV-*A +%7rG f 7rAr(-E')r 7r JV-»A(-E')) (4) 

Ta^ttn(E) = [1 + t 7TN (E)G 7rN (E)]T A _ r7rN . (5) 
The A self-energy in Eq. (J2J) is then calculated from 

^a(E) = T 7rN _ rA G 7rN (E)T A ^ wN (E). (6) 

As seen in the above equations, an important feature of the dynamical model is that the 
bare vertex T 7 7v^a is modified by the non-resonant interaction t> 77r to give the dressed vertex 
f 7 ./v— >a , as defined by Eq. (j3J). The second terms in the right-hand sides of Eqs. (4) and 
(5) can be interpreted as the dynamical pion cloud effects on the N-A transitions. We can 
then identify the parameters of the bare vertex T-^^a with the predictions from hadron 
models within which the 7r-baryon states in continuum are excluded. Such a separation of 
the reaction mechanisms from the excitations of hadron internal structure is an essential 
ingredient of a dynamical approach, but is not an objective of the approaches based on 
either dispersion relations fl2l k| or K-matrix method 14, 

In the above formulation, the matrix elements t> 77r and r 7 jv+-»A are determined by the 
electromagnetic current JJf 1 . The dynamical formulation for investigating weak pion pro- 
duction reactions can be obtained from Eqs. (l)-(6) by replacing J^ m by — A^, where 

and are the weak vector and axial vector currents respectively. By Conserved Vector 
Current (CVC) relations, the vector currents can be obtained from the electromagnetic 
currents by isospin rotations. Guided by the effective chiral Lagrangian method and using 
the similar procedures of the SL model, we can construct the axial vector currents A^ for 
pion production reactions. The resulting axial current amplitude consists of nucleon-Born, 
rho-exchange, pion-pole, and A excitation terms and can be calculated using the parameters 
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in the literatures. In particular, the bare axial N-A transitions at Q 2 — > are assumed to 
be those calculated by Hemmert, Holstein, and Mukhopadhyay[19] within the constituent 
quark model. 

Most of the earlier theoretical investigations 0, HE HE EE IHl of weak pion production 
reactions were performed during the years around 1970. The model which was most often 
used in analyzing the data was developed by Adler[25|. He considered a model based on 
the dispersion relations of Chew, Goldberger, Low, and Nambu(CGLN) The driving 

terms and subtraction terms of the dispersion relations are calculated from the nucleon Born 
terms. No vector meson exchanges are included. Additional subtraction terms are added 
to remove the kinematic singularities of multipole amplitudes. In solving the dispersion 
relation equations, only the P33 phases are included to account for the unitarity condition 
via Watson Theorem and the imaginary parts of other multipole amplitudes are neglected. 
Clearly, Adler's approach needs improvements for a precise extraction of the axial N-A from 
factors from the data. We however will not address this issue in this paper. Instead, we 
focus on the development of a dynamical approach outlined above. 

In section II, we present a formulation for calculating the cross sections of neutrino 
induced pion production reactions. The current matrix elements needed for our calculations 
are given in section III. The results are presented in section IV. Conclusions and discussions 
on possible future developments will be given in section V. 



II. CROSS SECTION FORMULA 

The formula for calculating the cross sections of neutrino induced pion production re- 
actions have been given in Ref.j25[. In this section, we would like to present a different 
formulation which is more closely related to the commonly used formulation of pion electro- 
production, and is more convenient for our calculations based on a dynamical formulation 
of the problem. 

We focus on the u+N — > l + n + N reaction, where I stands for electron(e) or muon(/i) and 
v for electron neutrino(z/ e ) or muon neutrino We start with the interaction Lagrangian 

L int (x) = Gf ^ S V (s)W + h -c-], (7) 
where G F = 1.16637 x 1(T 5 (GeV)~ 2 , cos9 c = 0.974, 

= $i(x)-f(l- ls )il> v (x) ( 8 ) 

is the lepton current and 

J*(x) = V(x) - A"{x) (9) 

is the hadron current. The vector current (V M ) and axial vector current(A M ) will be con- 
structed later in terms of hadronic degrees of freedom. 

The coordinate system for calculating the differential cross sections of the v{p v ) + N{p) — > 
KPi) + tt(&) + N(p') reaction is chosen such that the leptons are on the x-z plane. The 
momentum transfer q = p v — pi, which is also the total momentum of the final tcN system 
in the laboratory system, is along the z-axis. The final ttN system forms a plane which has 
an angle cj) n with respect to the lepton v-l plane. This choice is identical to that widely 
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used in the study of pion electroproduction. It is then straightforward to show that the 
differential cross section can be written as 



da 5 



dEidnidn% 



Gp cos 2 e c 


Pi 




k c \m N 


2 


Pu 


(2tt 


fW 



L^W^, (10) 



where is the nucleon mass, and the phase space factor is expressed in terms of lepton 
momenta pi and p u in the laboratory(Lab) frame, pion momentum k c in the irN center 
of mass(C.M.) frame, and the invariant mass W of the irN subsystem. The angle f2* in 
Eq. (jl())l is the pion scattering angle defined in the ttN center of mass frame. 
The lepton tensor in Eq. (JTUJl is 

IT = fxVl + PM - gTipu ■ Pi) ± ie a ^p v , aPl ,p. (11) 

Here e 0123 = 1 and +(— ) in the last term is for neutrino (anti-neutrino) reactions. The 
hadron tensor is 

WW = \ E wl* ( 12 ) 

S ' N S N 

where implicitly represents the matrix element 

1 



< k,ps N I J M | ps N > 



m N 



(2tt) 3 \ 2E 7T (k)E N (p')E N (p) 



J, (13) 



with Sat denoting the z-component of nucleon spin. Such a simplified notation for hadron 
current matrix element will be used in the rest of this paper. 

To see the current non-conserving part of hadron current explicitly, we introduce the 
following variables 

K" = P »+ p i , (14) 

q»=p»-pl (15) 
Q 2 = -tf = -m ! l+2(p v - Pl ), (16) 

where g M is the lepton momentum transfer and mi is the mass of lepton /. We then obtain 

±ie a ^q a K^jl\. (17) 

The hadron current in Eq. p7|) is defined in the Lab frame. To account for the final 
ttN interaction within the dynamical formulation and to get an expression in terms of the 
scattering angle f2* of Eq. (jlOj) . we need to write j M in terms of the current j£ defined 
in the irN CM. frame. With the choice of the coordinate system described above, the 
transformation from the Lab frame to the ttN CM. frame is just a Lorentz boost along the 
z-axis. For the x and y components of the hadron current, we obviously have the following 
simple relations 

f = cos^-sin^, (18) 
f = sin^ + cos^. (19) 
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For the z and time components, the transformation relations are 



where 



j" = coshaj° + sinhaj^, (20) 
j z = sinhaj° + cosher/*, (21) 



coshe* =^ =^±^, (22) 

m N W 



M = \1 

Here we have introduced g M = (u, q) and = (u c , q c ) as the momentum transfers in the 



sinha = i3sL = m. (23) 



Lab and nN CM. frames respectively, and E NjC = yq 2 + m%- 

Since the leptons are on the x-z plane and the momentum transfer q is in the z-axis, we 
obviously have 

(24) 
(25) 
(26) 



K x = 


\Pu\\Pl\ ■ a 

— i — i — sin &i 
\Q\ 


K y = 


o, 


K z = 


\Pv\ 2 - \Pl\ 2 


2|g| ' 


K° = 


\Pu\ + Ei 
2 



(27) 

In the above equations, d% is the angle between the outgoing lepton momentum pi and the 
incident neutrino momentum p u . 

By using the above relations Eqs. ()18| ) -()27| ) . it is straightforward, although tedious, to 
cast Eq. (|17j) into the following form 

L^Wp, = - J21 r t + Rl + Rlt cos 0, + R TT cos 20, 

s' N ,s N 

+R LT : sin 0, + R TT , sin 20,] (28) 

All terms in the above equation can be calculated from jg defined in the 7rN CM. frame. 
Explicitly, and Rl terms are given as 

Rt = ((K*) 2 + ^^XUcT + \f c \ 2 ) =F xlm(f c jr), (29) 



1 

+(y 2 -u; 2 + m 2 )\p\ 2 ], (30) 



Rl = —i(x 2 -Q*-Q 2 mt)\j°\ 2 + 2(xy-Q 2 u c -uj c m*)Re(jy) 



with 



(Ey + EiW-um? (31) 
y _ {E v + Ej)u c + cosh a mf 
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In deriving Eq. (|5Uj) , we have defined p = q ■ j 
The other terms in Eq. (28) are given as 



UcJc 



\q c \jc an d eliminated 



Rlt 
Rtt 
Rlt> 
Rtt' 



2K a 



;-Re((x Jc ° + yp)jT) ± Im((Q 2 Jc ° + co^jf) 



{K x ) 2 {\f c \ 2 -\3l\\ 
1K X 

-— (Re((*j c ° + yp)jT) ± Im((Q 2 Jc ° + u; cP )f c *), 



-2{K*fRe{f c f c *). 

It is interesting to note that in the mi = limit the above formula become 

Q 2 M 2 + W 



Rt 
Rl 

Rlt 
Rtt 
Rlt 1 
Rtt 1 



l-e L 2 

Q 2 Q 2 r 



TVT^ImO-if)], 



0|2 



1 - e 19c 

Q 2 



2 > Jc ' ' 



1 



2e(l + e)Q 2 



[-Re(j° c f c *) ± 



1 + e 



Im(J c iD], 



Q A\j x c \ 2 -\j y c 



1 - e 1 

q 2 

1 - e 

12 



\ 



<3 



1 -e 1 



2e(l + e)Q 2 



-eRe(^T)], 



[Re(J c °jD ± 



-MfJT)\ 



where 



1 + M- tan 2 



(33) 
(34) 
(35) 
(36) 

(37) 
(38) 

(39) 
(40) 
(41) 
(42) 

(43) 



with j° = j° + u c p/Q 2 . Eqs. (37)-(43) are very similar to the familiar forms of pion electro- 
production. In the energy region we are considering, the muon mass cannot be neglected. 
Calculations of the differential cross section Eq. ()10|) for the v^ + N— > p + n + N reactions 
must be done by using Eqs. (|29 p -()36j ) . 

To deal with the old data, we need to calculate the differential cross section in terms of 
the variables ^(invariant mass of ttN subsystem), Q 2 = — g 2 (Eq. ([16))) instead of lepton 
energy (Ei) and scattering angle Q\. The relation between them are 



Q 2 = 2E u E l - 2\p l \\p„\ cos -m 2 , 



W = yjm% + 2m N {E v - E t ) - Q 2 . 
By using the above relations, we can obtain the following relation 



da 



2irW 



da 



dWdQ 2 dtt* 7r 2m N \p u \\pi\ dEidHidn* 



(44) 
(45) 



(46) 



7 



To calculate the differential cross sections in the right-hand-side of the above equation using 
Eq. (|1U|) and Eqs. (|29 p -([36 p . we also need to know the variables in nN CM. frame. These 
are given by 

a, = W2 -t~ m \ (47) 



2W 

W 



\Qc\ = -firkl, (48) 



\^ = i( W2 + ^~ m2N f-<. (49) 
The total cross section for a given incident neutrino energy E v is then calculated by 

fWmax rQmax „ d(J 

a(E v ) = dW d Q 2 -^—, 50 

Jw min JQ 2 mm dWdQ 2 



where 



dWdQ 2 J n dWdQ 2 dni 



G 2 F cos 2 9 c 1 \k r \ r 1 



dni- £ [Rt + Rl]. (5i) 



2 327r 4 |p v | 2 7 ^2 

The integration ranges in Eq. ()50|) are found to be 

W min = m N + m n , (52) 
W max = W T - m h (53) 
where m n , m^, and mi are the masses of the pion, nucleon and the outgoing lepton respec- 
tively, and Wt = sj (Pu + Pn) 2 = ^2m^E v + m 2 N is the invariant mass of the initial v-N 
system. For a given allowed W, the range of Q 2 is found to be 

Qmin = -m 2 + 2E;(E?-ti), (54) 
QL X = -m 2 + 2E:(E? + p c), (55) 



where Ef = ymf + pf 2 , E c v and p\ are the neutrino energy and outgoing lepton momentum 
in the CM. frame of the whole system(not the CM. frame of the final ttN subsystem). 
Explicitly, we find 

K - « («) 

(58) 



III. CURRENT MATRIX ELEMENTS 



To proceed, we need to construct current operators = V^ — A^. The matrix elements of 
these currents are the input to solving the dynamical equations that are of the form of Eq. (2)- 
(6) with appropriate changes of notations; namely v^n — > vvn or van and Ta^ 7 v — > ^a->vn 
or Tj\_,an- In the first part of this section, we present the non-A current matrix elements. 
The amplitudes associated with the A excitation will be given in the second subsection. 
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A. Non-A Amplitudes 



We first consider the hadronic currents associated with ic and N degrees of freedom. Here 
the chiral symmetry is the guiding principle. These currents can be derived by using the 
well-developed procedures^, 3(|. Within the SL model, we need to also consider currents 



associated with p and lo mesons. Since the electromagnetic current(J^ m ) is related to the 
isovector vector currentfV^) by J^ m = j™°- scalar _|_ y% W e can use the Conserved Vector 
Current (CVC) hypothesis to obtain the weak vector current from Vj% of SL model by isospin 
rotation. We find 

2mjv 2 m n 



+ [tt x ^tt] ■ v, + ^e a ^[d Q A ■ *[8nu% (59) 



where is an arbitrary isovector function. Note that p meson does not contribute to the 
charged vector currents considered in this work since p-ir current is isoscalar. 

To construct axial vector currents associated with 7r and N degrees of freedom, we are 
guided by the standard effective chiral Lagrangian methods |2flL I30I] and follow the procedure 
of SL model. We then obtain the following form of axial vector current 

& = QaN^J-N - f pnA fT x 7? - Fd^ir. (60) 

Here F = 93 MeV is the pion decay constant, and g A = 1-26 is the nucleon axial cou- 
pling constant. Note that uj-tt current is G-parity violating second class current and is not 
considered here. 

With the above current operators and the following Lagrangians from SL model 

Ataw = -—NwsfN.d**, (61) 

L pN „ = g P W - i^°v>V?\ -^N + Tcxd^-fr), (62) 

Lujnn = 9u,nnN[0 - a^P^N, (63) 

zmjv 

we can evaluate the current matrix elements of the tree-diagrams illustrated in Figs. 1-2. 
For vector current contributions(Fig. 1), we have 



m 



V»{k,j,q,i) = ^ ^ 2K[k)EN{p)EN{pl) U{P') [VLrnik, 3, g, i) + V&(k, j, q, i))u(p) , 

(64) 

where the Born term(Figs. l(a)-(d)) is 

VBorn(k,j,q,i) = i[T 7T (k,j)S F (p + k)Vj}(q,i) + Vfi(q, i)3 F (p - k)T„(k,j)] 

+ ^ tT ,_ 7V __(^ ( , +p _ pr] (65) 
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with 



K(kJ) = — (66) 
V£(q,i) = [V + i^^ly. (67) 



and 



^ — rriN 

is the Dirac propagator. The u-n term (Fig. le) is 

y, m = + ^y-rt /^-^ . (68) 

The p — 7r term has a similar form, but it is a isoscalar and does not contribute the charged 
currents considered in this investigation. 

The non-A axial vector current amplitude has two parts. The first part is due to the 
first two terms of Eq. (60) and is illustrated in Fig. 2. The second part is due to the 
term —Fd^n. Obviously, this interaction will induce a pion-pole term illustrated in Fig. 3, 
where the shaded box is identical to Fig. 2 except that the axial vector field(waved line) 
is replaced by the pion field(dashed line). For the non-pion pole part of the axial current 
contributions(Fig. 2), we have 



(69) 

where the Born term is 

A" Born (k,j,q,i) = i[r n (kJ)S F (p' + k)A^ N (q,i)+A%(q,i)S F (p-k)r n (kJ)}, (70) 

with 

A%(q,i) = 9a1^ T -. (71) 

The p-ir term in Eq. (69) is 

A^AKj^i) = -9 P f P „ A e ijk ^[Y + V - PU Z _ p)2 ■ (72) 

Taking a phenomenological point of view, we fix the coupling constant f fm A in the above 
equation by using the soft pion limit 

in? 

fpnA = ■ (73) 

* 9o 
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The pion pole term(Fig. 3) can be easily obtained by modifying A^ P , defined by Eq. (69), 
to include a pion propagator. By using the PCAC relation, we find that this pion pole term 
can be easily included by the following procedure 

Aufi ■ -\ & (i ■ -\ q^q- A N p{k,j,q,i) , v 

Ar[k,],q,i) = A% P (k,j,q,z) = = . (74) 

qz _ m z 

In the dynamical approach, as briefly reviewed in section I, the non-resonant amplitudes 
will be integrated over the nN scattering amplitudes. Thus the non-A amplitudes, as il- 
lustrated in Figs. 1-2, must be regularized by introducing a form factor at each vertex. 
Furthermore, the finite size effects of hadron structure also require including form factors. 
Fortunately, these form factors can be taken from the SL model and other theoretical in- 
vestigations. The vector current matrix elements are regularized as those in SL model. For 
axial current matrix elements, the A-NN vertex is regularized by a dipole form factor 

Fd(q2) = + (75) 



where Ma = 1.02 GeV is taken from Ref. [31|. The n-NN and p-NN form factors are taken 
from SL model^. The A-p-ir form factor is also assumed to be of the dipole form of Eq. (|75j). 
With these specifications, the non-A amplitudes do not have any adjustable parameters in 
our investigations. 



B. A amplitudes 

The A has two contributions to weak pion production reactions, as shown in Fig. 4. 
The resonant amplitude is due to the formation of a A in s-channel(Fig. 4a). The cross 
A term, Fig. 4b, is part of the non-resonant amplitude. Each term has a corresponding 
pion pole term illustrated in Fig. 3 and these pion pole terms must be also included in our 
investigations. 

To calculate these A amplitudes, we need to define the matrix elements < A | \ N > 
and < A | A^ | N >. The vector current matrix element < A | V 1 \ N > can be obtained 
from SL model by appropriate isospin rotations. Here we focus on the axial vector current 
matrix element. 

It is well known 3, 25| that the most general form of axial vector current matrix element 
can be written as 



< A\A^\N >= u^( P ')Y^Tu{p) (76) 



with 



where u& u (p) is the spin 3/2 Rarita-Schwinger spinor, q = p' — p, P = p' + p, and T % is the 
i— th component of the isospin transition operator(defined by the reduced matrix element 
< 3/2 || T || 1/2 >= - < 1/2 || T+ || 3/2 >= 2 in Edmonds convention 0). 
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It is useful to explore here the meaning of each form factor of Eq. (|76|). Since u Av {p') 
u (p')tA,v(p'), where €a,u(p') is an unit four vector, we can rewrite Eq. (fTtjj) as 



< A\A"\N > = u(j>')[(dl(<?) + ^-d 2 (q 2 ))e A + (-d 2 (q 2 )P» + d,{q 2 )q^ ' 



m" N m* N 
+i^p-e^e A>v P a qp l5 ]T i u(p). (77) 



m 



N 



In the rest frame of a A on the resonance energy( p' = (toa,0),p = (E^(q),—q),q = 
(m A - E N (q),q),e A = (0, e), u(p' = (m A ,0)) = (x£,0), and u (p) ^ Xn, where \p is the 
Pauli spinor), the space and time components of Eq. (77) become 



< W > = + =^4)S - (4 + 4)^5 

V AiilN N rn N 

-id A s * q{(T - q \ ]r, (78) 

m 2 N (E N + m N y 

< A\A"\N > = lE N + mN [d2 S- q (mt + E N ) _ • q{m* - E N ) ]T ^ (?g) 

V 2?77,jv ^jV m ./V 

where we have defined the transition spin S = Xa € Xn (the reduced matrix element < § | 
S || g > is identical to that of the transition isospin T). The above expression suggests 
that di,d 2 terms describe the Gamow- Teller transition and d^ describes the quadrupole 
transition. The term of d^ comes from the pion pole term illustrated in Fig. 3. 

For simplicity, we now follow Ref.^J to fix the form factors di(q 2 ) at q 2 = using the 
non-relativistic constituent quark model. The axial vector current operator for a constituent 
quark is derived from taking the non-relativistic limit of the standard form gAqQl^l^Q- By 
using a procedure similar to Eq. (74) based on the PCAC relation, we then can extend the 
resulting current operator to also include the pion pole term contribution(Fig. 3) induced 
by the current —Fd^ir of Eq. (60). The total axial vector current operators associated with 
a constituent quark model are found to be of the following forms 

<p'\A l \p> = g Aq [a + Q a ■ q}^, (80) 

qz _ m z 2 

2m q q M — m£ 2 

where p is the quark momentum, o and r are the quark spin and isospin operators respec- 
tively. By using the standard non-relativistic s-wave quark wave functions for N and A, it 
is straightforward to obtain 

< p' N | A* | p N > = {^g Aq )[o- + q2 q m2 <r ■ Q] \, (82) 
<Pn\A'"\p k > = (§ ,M)£^± + ^ • 4 < 83 > 



and 




7 V HC-l q 
25 ( 3^ )[5+ ^"^ 



<p A \A*\p N > = - x l — (- gAq )[S+ 2 ? 2 S-q]T\ (84) 
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In above equations, Pn and p& are the momenta for N and A respectively, and we have 
expressed the right hand sides in terms of usual nucleon spin and isospin operators and the 
A-iV transition operators S and T(as defined above). 

If we set 3m g = m N and \$Aq — 9a(®) — 1-26, Eqs. (82) and (83) can be identified with 
the usual non-relativistic form for the nucleon axial current matrix element at q 2 = 0. This 
then fixes the value of g^q which also determine Eqs. (84)- (85) for the axial N-A transitions. 

Comparing Eqs. (84)-(85) evaluated at the A rest frame(pA = 0) and Eqs. (78)-(79) 
taken in the non-relativistic limit (En — > tun and q 2 /m 2 N — > 0), we find that 

^)^^)d + 2m ^- A ' + t K) ) (86, 

*«® = (87) 

d 3 (Ql) = sKQD-^-i (88) 

where g\(Ql) = ^^9a- Here we also follow Ref.fl^ to assume that the quark model 
results Eqs. (84)-(85) are for the momentum transfer q = (m\ — m^, 0) and Ql = —q 2 = 
(rriA — rn N ) 2 . Eqs. (86)-(88) agree with the results of Ref. [l9| if we neglect the difference 
between and tua- 

To account for the g 2 -dependence, we assume that 

d t (Q 2 ) = d l (0)E(Q 2 ) (89) 

where F(Q 2 ) will be specified in the next section. For a given choice of F(Q 2 ), we can use 
Eqs. (86)-(88) to obtain ^(0) = d(Q 2 ) / F (Q 2 ) . The form factors di(Q 2 ) for i = 1,2,3 of 
the covariant form Eq.(76), which is used in our numerical calculations, are then completely 
fixed by the constituent quark model calculation. We neglect the A deformation in this 
work and set d^(q 2 ) = 0. 

With the axial N-A vertex specified above and the ttN model constructed in Ref. [lj, the 
dressed AN — > A vertices (i.e. Eqs. (4)-(5)) can be calculated for evaluating the resonant 
amplitude Fig. 4a(i.e. the second term of Eq. (2)). With the 7riVA interaction Lagrangian 
L n NA = f-KNA/ m ni J A(. x )TN(x) ■ d fl 7i(x) of SL model, the cross A term (Fig. 4b) is found to 
be 



IlL(kJ,q,i) = i^^T*rS V 6(p-k)T j k*. (91) 



with 



where the Rarita-Schwinger propagator S IMU (p) given explicitly in Eq. (3.18) of Ref. |]|. 

The above derivations allow us to use the constituent quark model to calculate the A 
amplitudes. This completes our derivations of the weak current matrix elements which are 
the input to the dynamical equations Eqs. (2)-(4) with the sub-index 7 replaced by either 
V or A. 
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IV. RESULTS AND DISCUSSIONS 



As presented in subsection III. A, the non-A amplitudes, illustrated in Figs. 1-2, do not 
have any adjustable parameters in this investigation. While more theoretical investigations 
may be needed to refine these non-resonant amplitudes, it may not be necessary at this time 
since the weak pion production data are still limited. Instead, we focus on the investigation 
of the axial vector N-A transition form factors defined in subsection III.B. 

To proceed, we first recall the N-A form factors introduced in SL model. For the bare 
magnetic Ml N-A transition, it was taken as 

G M (Q 2 ) = G M (0)R SL (Q 2 )G D (Q 2 ), (92) 

where G D (Q 2 ) = 1/(1 + Q 2 /M 2 ) 2 with M 2 = 0.71 GeV 2 is the usual dipole form factor of 
the nucleon, and the correction factor is defined as 

R SL = (1 + aQ 2 ) exp(-bQ 2 ). (93) 

By fitting the pion photoproduction and electroproduction data up to Q 2 = 4 (GeV/c) 2 , 
it was found that G M (0) = 1.85, a = 0.154 (GeV/c)~ 2 and b = 0.166 (GeV/c)~ 2 . Here 
we remark that if we follow the same procedure given in section III.B to also calculate the 
7iV — > A transition using non-relativistic quark model, we find that the quark model yields 
G*m (Qo) = 2-37 for Ql = (m& — mN) 2 - This value is very close to the value Gm{Q1) = 2.42 
which can be obtained from SL model by using Eqs. (92)-(93). This suggests that the use 
of Eqs. (86)-(88) for defining di(Ql) of the axial N-A form factors is a good starting point 
of our investigations. As a continuation of SL model, we therefore first consider a model, 
called Model I, which assumes that the form factor of Eq. (89) is 

F{Q 2 ) = R sl (Q 2 )Ga(Q 2 ), (94) 

where R SL (Q 2 ) is given in Eqj93), and G A (Q 2 ) = 1/(1 + Q 2 /M 2 A ) 2 with M A = 1.02GeV of 
the nucleon axial form factor [3 1|. 

We first compare the total cross sections predicted by Model I with the dataj33[. The 
results are shown in Fig. 5. We see that the predictions(solid curves) agree reasonably well 
with the data for three pion channels. For the data on neutron target, our predictions(solid 
curves in the middle and lower figures) are in general lower than the data. This is perhaps 
related to the procedures used in Ref.|33j to extract these data from the experiments on 
deuteron target. 

One of the main features of the dynamical approach taken in this work is to renormal- 
ize(dress) the N-A transitions with the dynamical pion cloud effect, as described Eqs. (4)-(5). 
The importance of this effect is shown in Fig. 5. We see that our full calculations(solid curves) 
are reduced significantly to dotted curves if we turn off the dynamical pion cloud effects. 
If we further turn off the contributions from bare N-A transitions, we obtain the dashed 
curves which correspond to the contributions from the non-resonant amplitudes(Figs. 1-3). 
Clearly, the nonresonant amplitudes are weaker, but are also essential in getting the good 
agreement with the data since they can interfere with the resonant amplitudes. 

We next compare the Q 2 -dependence of the differential cross sections da/dQ 2 predicted 
by Model I with the data from ANL@ , BNL|2l|, and CERNQ. Here we need to account 
for the variation of neutrino flux in the experiments at ANL|33|| and BNLjijJ. We calculate 
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the following quantity 



a _ , f*~ ^ N(K) f- ^ N(E) 



dQ 2 J Emin @ model 

(E v ) dQ 2 v " n ' 

J Emin O model 

where N(E U ) is the distribution of events in neutrino energy E v which is within the range be- 
tween E min and E max , and a m nH.p.i(E„) is the calculated total cross section. The distributions 
N(E U ) are given in Fig. 6 of RefQ and Fig. 4 of Ref.Q. 

The predictions from Model I are compared with the ANL data|33| in Fig. 6. We see 
that our predictions(solid curve) agree reasonable well with the data both in magnitude and 
Q 2 — dependence. In Fig. 6 we also compare the contributions from vector current (dot-dashed 
curve) and axial vector current (dotted curve). They have rather different Q 2 -dependence 
in the low Q 2 region and interfere constructively with each other to yield the solid curve of 
the full results. Since vector current contributions are very much constrained by the (e, e'n) 
data, the results of Fig. 6 suggest that the axial vector currents constructed in section III 
are consistent with the data. 

For comparing with the BNL datajijj, we normalize the calculated da /dQ 2 to the events 



data of Fig. 5 of Ref. [21] at Q 2 = 0.2 (GeV/c) 2 . As shown in Fig. 7, the data can be 
reproduced very well by Model I. The BNL data was used in the most recent attempt |2l[ 
to extract the axial N-A form factor. We will discuss this later. 

The comparison with CERN dataj34| is given in Fig. 8. Here the data have some structure, 
which is perhaps mainly due to the poor statistics of experiment. We see that model I can 
reproduce the main feature of the Q 2 - dependence. 

To explore the effects due to the dynamical pion cloud on the N-A form factors, it is 
instructive to recall here the results of Ref.[2j for the magnetic Ml 7iV — > A transition. This 
result is displayed in the left panel of Fig. 9. We see that the pion cloud effect is essential 
in explaining the empirical values extracted directly from the data. The corresponding 
dynamical pion cloud effect on the axial N-A form factor is shown in the right panel of 
Fig. 9. We again see fairly sizable contribution from dynamical pion cloud. We stress that 
the empirical form factor, such as that extracted|2l| from BNL data, can only be compared 
with our dressed form factor, since the dressed form factor, not the bare form factor, directly 
determines the reaction amplitude Eq. (2). The differences between the solid and dotted 
curves in Fig. 9 explain the observation [19] that the quark model prediction of axial vector 
N-A strength at Q 2 — > is lower than the empirical value of Ref. [2j| by about 35 %. This 
is simply due to the fact that the dynamical pion cloud is not included in constituent quark 
model calculation and most of the hadron model calculations of N-A transitions. It is also 
interesting to observe from Fig. 9 that the pion cloud effect on the axial N-A form factor 
is mainly to increase the magnitude, not much to change the slope. On the other hand, 
both the magnitude and slope of the magnetic Ml form factor (left-hand side of Fig. 9) are 
significantly changed by including the pion cloud effects. 

To see the dependence of our predictions on the parameterization of axial N-A form 
factor, we next consider Model II which differs from Model I only in replacing the correction 
factor R sr, o f Eq. (94) by a different form. Rather arbitrarily, we consider a form used in 
Refs. [Til 21, 25| . The axial N-A form factor in this model II is also defined by Eq. (89) with 

F{Q 2 ) = R H {Q 2 )G A {Q 2 ) (96) 

where Ga{Q 2 ) is the axial nucleon form factor of Eq. (94) and 

R » = ( 1 + a TTQ^ < 97 > 
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with a = -1.21 and b = 2 (GeV/c) 2 . 

The main difference between Model I and Model II is in determining ^(0) = d(Ql)/F(Ql) 
from the quark model values di(Ql) given in Eqs. (86)-(88). We find that l/F(Ql) = 1.20 
for Model I and 1/F(Qq) = 1.26 for Model II. Thus these two models are already different 
at Q 2 — > 0. Their Q 2 -dependence are compared in Fig. 10. We see that the form factor of 
Model II(dotted curve) are much lower than that of Model I(solid curve). Consequently, the 
weak pion production cross sections calculated from Model II are found to be about 20 % 
lower than the data. This is illustrated in Fig. 11. We find that the model II can fit the 
data better if we increase the strengths di(0) of axial N-A form factor by 20%. The results 
shown in Figs. 10 and 11 demonstrate the sensitive of the data to the form of the bare axial 
N-A form factor. The data clearly favor the form factor Eq. (94). 

As seen in Eqs. (2) and (4), the empirical N-A form factor extracted directly from the 
data, such as those obtained in Ref . j2l| can only be compared with our dressed form factor 
F. The comparison is given in Fig. 12. We see that they do not agree well, in particular 
in the high Q 2 region. The differences mainly come from the fact that the non-resonant 



amplitudes of the Adler Model, which was usedj21| in extracting the empirical N-A form 
factors, are rather different from what we have in our dynamical model. With the very 
limited data, it is not possible to exclude one of these two rather different results. Clearly, 
more precise data are needed for a better test of our models. 

To stimulate future experimental efforts, we now present two predictions of Model I in 
Figs. 13 and 14. In Fig. 13, we present the Q 2 -dependence of differential cross sections 
at W = 1.1,1.2,1.236,1.3 GeV. The ^-dependence at Q 2 = 0.1,0.5,1.0,1.5 (GeV/c) 2 is 
shown in Fig. 14. 



V. CONCLUSIONS 

We have extended the dynamical model of Refs. 0, to investigate weak pion produc- 
tion reactions in the A region. The calculations for neutrino-induced reactions have been 
performed by using the current matrix elements illustrated in Figs. 1-3. The parameters for 
the non-resonant terms(Figs. 1-2) are taken from the previous investigations of pion electro- 
production. The bare N-A transitions at Q 2 — > are fixed by the constituent quark model. 
With the axial N-A form factor given in Eqs. (86)-(88), the predicted total cross sections 
and the differential cross sections da/dQ 2 agree reasonably well with the existing data. 

We have analyzed the calculated N-A form factors. Similar to the finding of the (e, e'ir) 
studies^, 0|, we show that the constituent quark model prediction of axial vector N-A form 
factor at Q 2 — ► is consistent with the data. The discrepancy observed in Ref. [l^] is due 
to the dynamical pion cloud effect which is not included in the constituent quark model 
calculation. 

An unsatisfactory aspect of the present investigation is the lack of extensive and good 



quality data. Thus the axial vector N-A from factor extracted [2Jj from the data relied 
heavily on model assumptions and hence the origins of the differences seen in Fig. 12 are not 
clear. It will be very useful to have sufficient data for performing partial wave decomposition, 
like what have been routinely performed in (e, e'ir) studies, such that the N-A form factor 
can be extracted model independently. Hopefully, the situation will be improved in the near 
future when new neutrino facilities will become available. 
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(a) (b) (c) (d) (e) 

FIG. 1: Non-A vector current contributions. Dashed lines represent pion. Waved lines represent 
the vector field. 



♦ » 




(a) (b) (c) 

FIG. 2: Non-A axial vector current contributions. Dashed lines represent pion. Waved lines 
represent axial vector field. 
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FIG. 3: Pion-pole term of axial vector current contributions. The dashed-box represents the either 
the non-A mechanisms of Fig. 2 or the crossed A mechanism of Fig. 4(b). 




(a) (b) 

FIG. 4: The A current contributions. Dashed lines represent pion. Waved lines represent either 
the vector field or axial vector field. 
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FIG. 5: Total cross sections of N{u^, (j,~ir)N reactions predicted by Model I. The data are from 
Ref . [33| . The solid curves are from full calculations. The dotted curves are from turing off pion 
cloud effects on N-A transitions. The dashed curves are the contributions from the non-resonant 
amplitude. 
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FIG. 6: Differential cross sections da/dQ 2 of p(z^ /i , /j,~Tt + )p reaction averaged over neutrino energies 
0.5 GeV < E u < 6 GeV. The curves are from calculations using Eq. (95). The dotted curve(dot- 
dashed curve) is the contribution from axial vector current A (vector curent V). The solid curve is 
from our full calculations with V-A current. The data are from Ref. 33 1. 
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FIG. 7: The differential cross sections da/dQ 2 of ,vr + ) reaction calculated using Eq. (95). 

The data are from Ref . 21 ] . 
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FIG. 9: The N-A form factors predicted by Model I: left panel: Magnetic Ml form factors given 
in Ref. 2], right panel: axial vector form factor. The solid curves are from full calculations. The 
dotted curves are obtained from turing off the pion cloud effects. Gd = 1/(1 + Q 2 /My) 2 with 
My = 0.84 GeV is the usual proton form factor and G A = 1/(1 + Q 2 /M\) 2 with M A = 1.02 GeV 
is the axial nucleon form factor of Ref. [3l|| . 
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FIG. 10: Bare axial N-A form factors for Model I(soild curve) and Model II(dotted curve). 
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FIG. 11: The total cross sections of fi ir + )p calculated from Model I(soild curves) and Model 
II (dot curves). 
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FIG. 12: Compare the dressed axial N-A form factor predicted by Model I (solid curve) with the 
empirical form factor (dot-dash curve) determined in Ref . [2l| . 
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FIG. 13: The predicted differential cross sections da/dWdQ 2 for p(z^, pT ir + )p at neutrino energy 
E u = 2 GeV and invariant mass W= 1.1 (dashed), 1.2(dotted), 1.236(solid), 1.3(dot-dashed) GeV. 
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FIG. 14: The predicted differential cross sections da/dWdQ 2 for p(z^, fi~n + )p at neutrino energy 
E v = 2 GeV and Q 2 =0.1(solid), 0.5(dashed), l.O(dotted), 1.5(dot-dashed) (GeV/c) 2 . 
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